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Abstract 

We investigate the properties of the single A hypernuclei within a Skyrme-Hartree- 
Fock (SHF) model. The parameters of the Skyrme type effective lambda-nucleon (AJV) 
interaction are obtained by fitting to the experimental A binding energies of hypernuclei 
with masses spanning a wide range of the periodic table. Alternative parameter sets are 
also obtained by omitting nuclei below mass number 16 from the fitting procedure. 
The SHF calculations are performed for the binding energies of the A single-particle 
states over a full mass range using the best fit parameter sets obtained in these fitting 
procedures and the results are compared with the available experimental data. The 
data show some sensitivity to the parameter sets obtained with or without including the 
nuclei below mass 16. The radii of the A orbits in the hypernuclear ground states and 
the A effective mass in nuclear matter show some dependence on different parameter 
sets. We present results for the total binding energy per baryon of the hypernuclei 
over a large mass region to elucidate their stability as a function of the baryon number. 
We have also employed the our best fit AN parameter sets to investigate the role of 
hyperons in some key properties of neutron stars. 

Keywords: Skyrme Hatree-Fock Model, Structure of Hypernuclei, Hyperon-nucleon 

effective interaction 
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1. Introduction 



Hypernuclei provide an unique opportunity to investigate the dynamics of the full 
meson and baryon SU(3) flavor octet. They are an excellent tool to extract informa- 
tion on the hyperon-nucleon interaction. For a complete understanding of the baryon- 
baryon interactions in terms of the meson-exchange or quark-gluon pictures, the basic 
experimental data are required also on the hyperon-nucleon interaction together with 
the nucleon-nucleon interaction. Yet, relatively little is understood about nuclei that 
contain one or more hyperons. It is possible that the study of hypernuclei may help 
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in unraveling fundamental issues regarding the ordinary nuclei in terms of quantum 
chromodynamics (QCD) which is the theory of strongly interacting particles. 

Since A is the lightest among the hyperons (toa - 1.115 GeV), the A hypernuclei 
are the most investigated systems sofar|[ll|3,[3lB|5lSI3,[l,l3- Through a series of 
experimental studies involving {K^,7t^), (tt^, K^), (y, K^) and (e, e'K^) reactions (see, 
e.g. Ref. 101 for a recent review of the experimental scenario) the A shell structure 
has been mapped out over a wide range of the periodic table. Systematic studies of 
the energy levels of light A hypernuclei have enabled the extraction of considerable 
amount of details about the AN interaction. It is established that the spin-orbit part 
of this force is weaker (see, e.g. Ref. |4] and the related references cite there) than 
that of the A^A^ system. This has also been useful in developing the theoretical models 
of the hypernuclear production reactions where the guantum numbers and the binding 
energies of A single particle states are vital inputs llOlll ill 12lll3LI 14111511 . 

For several reasons it is necessary to have information about the behavior of a 
hyperon in nuclear medium which can be provided by the heavier A hypernuclei. The 
basic quantities like A effective mass and the depth of the A potential in nuclei can be 
obtained by investigating the A binding energies in a wide range of heavier systems. 
These quantities are the important basic parameters for a realistic discussion of neutron 
stars 11611 . The information on the A spin-orbit splitting in heavier hypernuclei is of 
interest because it may have contributions from the higher order many body effects in 
addition to the AA^ two-body spin-dependent interaction. 

It is therefore, important to develop reliable theoretical tools to investigate the spec- 
troscopy of the hypernuclei. From QCD point of view hypernuclei lie in the non- 
perturbative low momentum regime. Therefore, lattice QCD calculations should be the 
ideal tool for investigating their structure. Indeed, first step in this direction has already 
been initiated where the scattering length and the effective range for AA^ scattering have 
been extracted in both QCD and partially-quenched QCD calculations iHtIi . However, 
a description of the detailed structure of hypernuclei is still beyond the reach of the 
lattice QCD and one has to use methods where baryons and mesons are the effective 
degrees of freedom. 

Both relativistic and non-relativistic descriptions have been used to investigate the 
structure of h ype rnuclei. The relativistic mean field approach has been employed in 
Refs. iflsjlTgiMIIllB with meson-hyperon vertices adjusted in varied ways. The 



SU(3) symmetric field theories including chirality 1123112411 . and the quark-meson cou- 
pling model LSi 25 , 3l that invokes the quark degrees of freedom have been developed 
to investigate the structure of hypernuclei. A density dependent relativistic hadron 
(DDRH) field theory was used in Ref. 15] to describe A hypernuclei. In Ref. [9] a in- 
medium chiral SU(3) dynamics has been used to study the hypernuclear single particle 
spectra where a very good agreement has been achieved with the corresponding data. 
The smallness of the A-nuclear spin-orbit interaction also finds a natural explanation 
within this approach. 

Among the non-relativistic approaches are the calculations based on the shell model 
? — I 

picture (see, e.g., Ref. 127] for a review). Most of them reproduce reasonably well the 
measured hypernuclear states of medium to heavy hypernuclei using phenomenolog- 
ical AN potentials of the Woods-Saxon type 0281 l29l ll 111 . Several authors have also 
developed the semi-empirical mass formulas to describe the binding energy (BE) of 



2 



the hypernuclei SIB in 

The Skyrme-Hartree-Fock (SHF) model, which has been a powerful tool for in- 
vestigating the properties of nonstrange nuclei ll33i [34l [ssll . was extended to studies 
of the A hypernuclei in Refs. |3^ 37]. The suitability of this approach for describing 
the A hypernuclei depends heavily on the proper knowledge of the AN effective inter- 
action, which has either been extracted from the microscopic methods Uke G-matrix 
calculations performed with Jiilich and Nijmegen potentials 14311 o r 

by fitting directly to the hypernuclear data (mainly the binding energies 
The reliability of the latter method is related directly to the amount of hypernuclear 
data used in the fitting process. In previous studies, the data used for this purpose were 
limited to a few light nuclei. Therefore, there is a need to perform SHF calculations 
using the AA^ interaction with parameters that are determined by fitting to a larger set 
of data on binding energies that is available now. 

In this paper, we study the A hypernuclei within the framework of the Hartree- 
Fock method where the parameters of the Skyrme type AN effective interaction have 
been determined by fitting to the experimental binding energies of of hypernuclei with 
baryon numbers ranging over a wide mass range. In one fitting procedure, we have 
included nuclei with masses ranging between 8 to 208. In two more fittings, the nuclei 
with masses below 16 were omitted. We have also included a density dependent term in 
the AA^ interaction in our fitting procedure as suggested in Ref. |,46i , ,47J . The parameter 
sets having the minimum values in three searches have been used to calculate the 
binding energies of about 95 hypernuclear states and the results are compared with the 
available experimental data. Furthermore, we have also calculated the total binding 
energy per baryon of 73 A hypernuclei, and the root mean square (RMS) radii of a few 
nuclei of them. Finally our best fit parameter sets have been used to investigate some 
properties of neutron stars. 

The present paper is organized as following. In section II, the Skyrme-Hartree-Fock 
method as applied to the description of the single A hypernuclei, is briefly described 
and parameters involved in the Skyrme type AA^ force are discussed. In section III, 
we describe the minimization method [based on the simulated annealing method 
(SAM)] that is used by us to determine the best fit parameters of the AA^ interaction. In 
section IV, we present the results and discussions of the comparison of our calculations 
for the single particle energies of the A states with the available experimental data. 
The results for the total binding energies per baryon of a larger number of single A 
hypernuclei are also presented in this section. We also investigate the role of our best 
fit AN interactions on some properties of neutron stars. Finally, in section V, we present 
the summary and conclusions of our work. 



2. Formalism 

The SHF model for nonstrange nuclei has been discussed in great details in Refs. fsl 



34113511 . In its extention 1 36, 37] to describe the A hypernuclei a contribution is added to 



the original energy density functional to account for the action of the hyperon-nucleon 
force. Thus, the total energy density functional (EDF) of a hypernucleus is written as 
a sum of two basic contributions - &n, which is the total energy density for neutron 
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and proton 113 311 and £a, which is the contribution due to the presence of the hyperon 
(hyperons): 



(1) 



where pg, Tq and Jq represent the baryon, kinetic and spin-orbit current densities, re- 
spectively (q = n,p,A). There are some additional terms contributing to the total hy- 
pernuclear energy that will be discussed latter on in this section. 

&N is related to the original SHF nuclear Hamiltonian density (Hn) as 



-I 



(2) 



where H^^ is written as 11481 \4' 



Hm — 



ft2 



1 



:^tn + Tf()[(2 + XQ)pij - (2x0 + l)(Pp + P„)] 



2m 

Y^t-.p^iil + x,)pI - (2x3 + 1)(pI +pI)] 

^[tl(2 + Xi) + t2(2 + X2)]tnPn + ^[f2(2x2 -I- 1) - ti(2xi + l)](TpPp + T„p„) 



+ ^[3fi(2 + xi) - t2(2 + X2)](Vp^)2 - ^[3fi(2xi + 1) + f2(2x2 + m(Vppf + (VpJ^] 
1 

+ 2 ^o[JiV-Vpiv + Jp-'^Pp + J«.Vp„] 



1 

16 



— (flXi + t2X2)ili + ^(fi - t2)[fp + il] 



1 2 CPpiyWr' 3 



4^ Ppir) 



'iPp(i-) 



1/3 



(3) 



where the first term on the right hand side represents the kinetic energy. We define, 
Pn - Pn+Pp, tn = T„+Tp and - Jn+Jp, where n(/?) correspond to a neutron (proton). 
We employ SLy4 Skyrme parameterization [48il in the calculation of the energy density 
of the hypernuclear core ^ 'Z. This interaction has been widely used in studies of the 



nuclear structure of normal and neutron rich nuclei OSOti . and the properties of nuclear 
matter and neutron stars |5l|] . 

The densities pq, Tq and are expressed as 



Tq(r) = I V(Pp(r,q) \\ 
3q{r) = ^ v^0*(r, ^)(-/V X cr)0^(r, q). 



(4) 
(5) 
(6) 



where 0/3(r, q) are the wave functions of the single particle states, vp represents the 
occupation probability, and cr the Pauli spin matrices. The sums are taken over all 
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occupied states for different particles q. Summations over spin, and isospin indices are 
implicit. 

The occupation probabilities vp of neutrons and protons are calculated by includ- 
ing the pairing energy functional in Eq. (1) as described in Ref. llsoll . The BCS 
equations for the pairing probabilities are obtained by variational method with respect 
to vp and written as. 
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1- 



(7) 



where ep is the single particle energy of the occupied state, and is the chemical 
potential. The pairing gap equation has the form - Gq Yj^eq ^jvp(_l - vp), where Gq 
is as defined in Ref. 15011 . 

The energy density functional £a is written as 



(8) 



As in Ref. BTI] . Ha in Eq.[8]has a term that corresponds to a Skyrme type two-body 
force. 



H 



AN 



= :^—ta + uq(1 + ^yo)PNpA 
ZniA 2 

+ ^(Mi + U2){Tf^pN + TnPa) + ^(3mi - "2)(VpN ■ Vpa) 
+ ^Wa(Vpn • Ja + VpA ■ Jn), 



(9) 



In ref. 113711 and II45I1 . a second term is added to Ha that corresponds to a zero range 
three-body ANN force. 



^3 2 

■jPa(Pn + '2-PnPp)- 



(10) 



This is completely analogous to the corresponding term of the nucleon Skyrme force 
proposed in Ref. f33\. It has been pointed out in Refs. |52, 53, Hij] that the binding 
energies of particularly light hypernuclei can be better reproduced if three-body AA^A^ 
term is included in the calculations. For example the overbinding problem of ^He 
hypernucleus has been solved by taking into account the AA^A^ force in Ref. 15511 . In 
the Skyrme Hartree-Foch calculations of Ref. ll45[l it is shown that the inclusion this 
term leads to a better reproduction of the Ip-ld and Is-lp levels spacings in several 
nuclei. 

In the Hartree-Fock calculations of even-even nuclei, the three-body part of the 
Skyrme interaction is equivalent to a two-body density-dependent interaction l33l 13511 . 
Therefore, some authors ll46ll47ll have included in Ha a term dependent on the nuclear 
density instead of the AA^A^ force. 



H 



ANp 



3 , 1 , 

-M3(l + -y3)p^ PA, 



y+1 



(11) 
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where y and are additional parameters. It may however, be remarked that at least for 
y + \, Eq. (11) can not be derived from from any three-body force. Furthermore, the 
equivalence of three-body and the density dependent forces in the hypernuclear case is 
not exact even for y - 1. In the limit of y = 1 and 3/3 = 0, Eq. (11) does reduce to 
a form similar (but not strictly equal) to Eq. (10). Nevertheless, using Eq. (11) leads 
to a value of nuclear incompressibility that is closer to its experimental value ll47ll . 
It is stated in Ref. |46] that Eq. (10) is not so adequate for representing the density 
dependence of the AA^ G matrix while Eq. (11) with a value of 7 = i is good enough 
for parameterizing the G matrix result. Therefore, in our calculations we have added 
Eq. (11) to the Hamiltonian instead of the three-body force term, Eq. (10), and have 
considered for y both 1 and |. We have used fp' jlnij^ - 17.44054 MeV fm^. 

Now, the single-particle wave functions 0;3(r, q) and corresponding single-particle 
energies ep are obtained by solving the Hartree-Fock equations with position and den- 
sity dependent mass term 



ft2 

-'^^TTTT-'^ + - ■ (V X cr) 

2m*{r) 



(r,q) = ep^p(r,q), (12) 



where m* is the effective baryon mass, and Vq and Wq are the central and spin-orbit 
terms of the mean field potential, respectively. The central term consists of (i) a purely 



nuclear part (Vnn) as described, e.g., in Refs. 11331 14811 . (ii) additional field created 
by the A hyperon (V^) that is seen by a nucleon, and (iii) the whole nuclear field 
experienced by the A hyperon (V^). The effective mass of the nucleon will also have 
additional terms. The spin-orbit part of the AA^ interaction has been ignored (Wa = 0) 



from the start as is done in Ref. Il37ll . The smallness of this term has been supported by 



several microscopic calculations iMEllal. The potentials and are written as 



= Mo(l + ^yQ)pN + ^(mi + M2)(Vta/5N + tn) 

+ i(3Mi - m2)(VVn) - ^(3mi - M2)(VpN/r) + ^M3(l + ^y3)pl^\ (13) 



and, 



= Mo(l + ;^3'o)Pa + 7(^1 + «2)(ta -H VtnPa) + ^(3mi - M2)(VVa) 

1 3 1 

- -(3mi - M2)(VpA/r) + -M^(l + -y^Xy + \)pIpk- (14) 

The spin-orbit part of the mean field has purely nucleonic contribution and it can be 
written in a straight forward standard way. 

The term containing the effective mass of the single A hyperon is expressed as, 

1 

TT^ = 7; + 7(mi + M2)PAr- (15) 

2m^ 2mA 4 

Similarly, the effective mass term for nucleon (in present of a hyperon) is written as 

tP- \ 1 1 

+ -[t\{2+xx)+t2(2+X2)]pN + -^[h(2x2+l)-t\{2xi + \)]pq, + -(ux+U2)pK, 



2mqr. 2mq. 8 ^ gL.v . . ,v , ,,r-q ^ 



(16) 
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where q' represents a nucleon (neutron or proton). The potentials, the effective masses, 
and the orbitals in Eq. ( fT2] l are evaluated alternatively until self-consistency is achieved. 

The total hypernuclear energy, in a density dependent Hartree-Fock model, includes 
contributions also from the H^nn and //^^ terms. Our total energy due to nucleon 
terms only includes relevant contributions. For the A case, we have 

fi« = J d\{\ + i3'3)pr PA (17) 

We also have to introduce to the total energy the center of mass (cm.) correc- 
tion arising due to the breaking of the translational invariance in the mean field of the 
Hartree-Fock theory. This is written as 

< p2 > 

(18) 



2(A — n)mM + nm\ 

where and mA are the masses of nucleon and A hyperon, respectively while n 
denotes the number of lambda particles (in our case n=l), and A is the baryon number 
of the hypernucleus. Pc.m. - 2* Pk is the total momentum operator in the cm. frame, 
which is the sum of the single particle momentum operators (pi^). Although the BCS 
state is not an eigenstate of Pc.m., and has vanishing total momentum < Pen. >- 0, it 
has non- vanishing expectation value of < P^.„ > given by, 

<Pc.,n.> = 2v^<aa|p2|aa> (19) 

a 

- ^ VttV/s(VaV/j - MaM;?) < ffyS I Pl.p2 I ay6 >, (20) 

a.l3 



where - - vl,, and a and yS represent the single particle states. This correction, 
however, is computed after variation (i.e., posteriori). The square of the single particle 
momentum operator appear only in the direct term of the cm. correction. The second 
and third terms represent the off diagonal single particle matrix elements of the mo- 
mentum operators that result from the exchange terms in < P^„, >. The cm. energy 
correction due to hyperons is approximated by the diagonal part of the cm. kinetic 
energy only. 

Thus, the total hypernuclear energy is given by the expression 

- £'NiPn,Pp, T„,Tp, J„,Jp) + SpairiVp, V„) + £f^ip„, Pp, Pf^, Ta) - £^ - (21) 

where Spai,-{vp, v„) is the pairing energy density functional as described above. 
The A binding energy in the SHF formalism is defined as 

Ba = ~ ^lA' (22) 

where £^ and £^ are the total binding energies of the hypernucleus and the core 
nucleus, respectively and binding energy per baryon number is obtained by dividing 
Ba with A (baryon number). 

The fitting procedure to obtain the parameters of the AN interaction is described in 
the next section. 
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Table 1 : Experimental values of the binding energies (in MeV) of various hypernuclei used in the x mini- 
mization procedure. 



Hypernuclei 




BE(Expt.) 






Is 


IP 


Id 


'He[59j 
iLi[59] 
]f Be [59] 


7.16 + 0.70 






8.50 + 0.12 






9.11+0.22 






]^"B [60, 61, 59J 


8.89 + 0.12 






aB[59,62J 


10.24 + 0.05 






i^B [59,63] 


11.37 + 0.06 






A C [64] 


10.76 + 0.19 






j^^C [65, 66] 


11.69 + 0.12 






i'N[67j 


13.76 + 0.16 






i^O [68,69] 


12.50 + 0.35 






a'SI [1 61J 


16.60 + 0.20 


7.0 + 0.2 




I'S [70] 


17.50 + 0.50 






^°Ca [71^72] 
i'V[73,Ml 


18.70+ 1.1 






19.90+1.0 






aY [7,6^ 


23.10 + 0.50 


16.50 + 4.1 


9.1 + 1.3 


TLa [7] 


24.50+ 1.20 


20.4 + 0.6 


14.3 + 0.6 


f Pb [7] 


26.30 + 0.80 


21.90 + 0.6 


16.8 + 0.7 



3. Parameterization for AN Skyrme Potential 

The values of parameters mq, ui, U2, u'^, yo and ^3 have been determined by fitting 
to the experimental binding energies (BE) of a set of hypernuclei across the periodic 
table by a minimization procedure that is based on the SAM. This is an elegant 
technique of searching for a global minimum in the hypersurface of the functions, 
which depend on the values of the parameters of the Skyrme interaction i58tl . This 
method has been used in Ref. [49,1 to determine the parameters of the Skyrme type 
nucleon-nucleon eff^ective interaction, which have been used to describe successfully 
the properties of normal and neutron rich nuclei as well as those of neutron stars. 

The;\f^ function is defined as 



Nd-N~ 



X 

V — 



1 



(23) 



where A^^ is the numbers of the experimental data points and A^^ the number of fitted 
parameters, cr, stands for the theoretical error and M*^^'' and M'!' are the experimental 
and the corresponding theoretical values, respectively, of a given observable (which in 
this case is the binding energy). The values of depend on the Skyrme parameters 
being searched [in our case mq, mi, M2, M3, yo and ^3 appearing in Eqs. (9) and (11)], 
since Mf is calculated using potentials with these parameters. Although there are 
correlations in some of these parameters (e.g. mq, yo, and m^, 3/3), we have varied 
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Table 2: The lower (vo) and upper (vi) limits, maximum displacement (d) and initial values (vjn) for the 
Skyrme force parameters used for implementing the SAM algorithm for minimizing the value that leads 
to the parameter set (HPA2) to be shown in Table[3]. 



Parameter 


Vo 


Vl 


d 


Vin 


Mo(MeV /m^) 


-592 


-275 


2.0 


-346 


Mi(MeV /m^) 





100 


1.0 


60 


M2(MeV /m^) 





100 


1.0 


50 


M;,(MeV fm^^^y) 


1850 


2300 


3.0 


1620 


yo 


-0.7 


-0.05 


0.01 


-0.13 


^3 


-1.5 


-0.01 


0.01 


-0.38 



Table 3: Parameters for the AN Skyrme potential derived self-consistently by fitting to the experimental 
binding energies of table I for two different values of y. 



SET 


r 


Uo 


Ul 


U2 


"3 

(MeV fm^+^r) 


yo 


y3 


x' 






(MeV fm^) 


(MeV fm^) 


(MeV fm^) 








HPAl 


1 


-326.395 


72.627 


-8.584 


1746.041 


-0.223 


-0.389 


3.58 


HPA2 


1 


-399.946 


83.426 


11.455 


2046.818 


-0.486 


-0.660 


3.06 


HPA3 


1/3 


-498.515 


65.203 


19.001 


995.832 


-0.492 


-0.444 


3.69 


HPA4 


1/3 


-475.584 


99.058 


-20.890 


1375.172 


-0.350 


-0.724 


3.76 



each one them independently - minimization of the x was the only requirement of 



our procedure. This is in complete analogy with the procedure of Refs. 115 811491] . The 
theoretical errors (see, Eq. l23l l in the fitting algorithm are taken to be equivalent to the 
experimental uncertainty in the data. 

The SAM is an elegant technique for optimization problems of large scale, in par- 
ticular, where a desired global extremum is hidden among many local extrema. This 
method has been found to be an extremely useful tool for a wide variety of min- 
imization problems of large non-linear systems in many different fields (e.g., see, 
Refs. Ill IzIItI)- Recently, the SAM was used to generate some initial trial pa- 



rameter sets for the point coupling variant of the relativistic mean field model 11771 17811 . 

In the SAM one needs to specify the appropriate annealing schedule together with 
the parameter space (i.e., the range of the values of the parameters) in which the best fit 
parameters are to be searched. As in Ref. Ii49il . we have employed a moderately faster 
Cauchy annealing schedule given by 

T(k) = TJk, (24) 

where T, is the initial value of the control parameter and T(k), with k= 1,2,3, , is 

the control parameter at the kth step. The value of k is increased by unity after l20Np 
reconfigurations or 12A^y, successful reconfigurations, whichever occurs first. The value 
of r, is taken to be 5.0, which is the same as that used in Ref. [49]. We keep on 
reducing the value of the control parameter using Eq. (l24l l in the subsequent steps until 
becomes stationary with further reduction. Since in our case number of parameters 
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Table 4: Experimental binding energy (in MeV) of various states of known single lambda hypemuclei with 
masses > 16 used in the alternative fitting procedure. 



Hypernuclei 






BE(Expt.) 








Is 


IP 


Id 


If 


Ig 




13.76 + 0.16 


2.84 + 0.16 










12.50 + 0.35 


1.85 + 0.05 










16.60 + 0.20 


7.0 + 0.2 










17.50 + 0.50 












18.70+1.1 










Sly 


19.97 + 0.13 


11.28 + 0.6 










23.10 + 0.50 


16.50 + 4.1 


9.1 + 1.3 


2.3 + 1.2 




24.50+ 1.20 


20.4 + 0.6 


14.3 + 0.6 


8.0 + 0.6 


1.5+0.6 


208pb 


26.30 + 0.80 


21.90 + 0.6 


16.8 + 0.7 


11.7 + 0.6 


6.6 + 0.6 



to be searched is rather large, we have defined the parameter space directly in terms 
of the range of values each parameter can take. Therefore, one of the key ingredients 
required for implementing SAM, in the present case, is to specify lower and upper 
limits for each of the parameter so that minimum is searched within these limits. In 
Table |2] we give some details of this procedure for one of the parameter sets to be 
shown in Table[3]that leads to the lowest In this table, lower and upper limits of the 
values of the Skyrme parameters are denoted by Vo and vi in each case, respectively. 
In the third column, d represents the maximum displacement allowed in a single step 
for a given Skyrme parameter during the reconfiguration process (see Refs. ll49ll and 
ifsill for more details). Vin, in the last column, shows the initial values of the Skyrme 
parameters used as a starting point for the SAM. 

We first applied this procedure to fit the BE data listed in Table [1] In our search 
procedure, we used two values for the parameter y, 1 and |. For each y, 2 parameter 
sets, generated by following the fitting procedure discussed above, are shown in Table 
3. The values for each set is given in the last column of this table. Sets HPAl and 
HPA2 have y = 1 while HPA3 and HPA4 have y - \- We note that parameter set 
HPA2 has the lowest j^^^ value among of all the sets. With y = 5, set HPA3 provides 
the lowest We have used the set HPA2 in further calculations of the hypernuclear 
properties unless specified otherwise. In neutron star studies, however, the set HPA3 
has also been used. The different sets of parameters shown in Table [3] lead to a depth 
of the AN mean field of about 28.00+0.58 MeV at the nuclear matter density (po) of 
0.16+O.Olfm-^ 

It could, however, be argued that hypernuclei with baryon number below 16 are too 
light to be used in the Skyrme Hartree-Fock fitting procedure, even through the Hartree- 
Fock method itself has been used to calculate the properties of light nuclei like He and 
Be by several authors (see, e.g. Ref. [93] and [41]). We have, therefore, also performed 
our fitting procedure by excluding all the nuclei with masses less than 16. Furthermore, 
since, experimental binding energies of If orbit in ^''Y and If and Ig orbits in ]^"'^La 
and ^'"^Pb are also available (see, e.g., Refs. [ 64. 7J), we have included them also in 
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Table 5: Parameters for the AN Skyrme potential derived self-consistently by fitting to the experimental 
binding energies of table|4]for two dift'erent values of y. The sets NAl and NA2 are obtained by including the 
BEs of 'j^N and those of the rest of the nuclei but excluding the BE of while OAl and 0A2 corresponds 
to the fits achieved by including BEs of '^O and the rest but excluding the BE of '^N 



SET 


r 


U{) 


Ul 


U2 


(MeV 


yo 


y3 








(MeV fm^) 


(MeV fm^) 


(MeV fm^) 








NAl 


1 


-253.3250 


147.1264 


-83.5843 


1684.9876 


0.5802 


0.4831 


2.30 


OAl 


1 


-236.5835 


116.8704 


-112.8812 


1453.3493 


0.1271 


-0.3110 


1.82 


NA2 


1/3 


-518.620 


82.0944 


-19.9772 


1190.1894 


-0.1392 


0.3126 


2.51 


OA2 


1/3 


-417.7593 


1.5460 


-3.2617 


1102.2221 


-0.3854 


-0.5645 


1.92 



our this alternative fitting. The BEs of the Is and Ip levels of the hypernucleus j^°N 
have recently become available from the study of the (e, e'K) reaction in the Jefferson 
Laboratory [67J while those for the j^^O hypernucleus are given in Ref. [ 73,1 . These 
binding energies together with those of the other heavier hypernuclear systems are 
shown in table H] 

We have performed two sets of fits - one by including BE data of the ]^''N hypernu- 
cleus together with those the rest of the nuclei shown in table|4]but excluding the '^O 
data, and the another where the BE data of '^O system are included but those of '^N 
are excluded. The resulting parameters having the minimum values, are shown in 
table|5]with gamma parameter of 1 and 1/3. The parameter sets NAl and NA2 include 
the binding energies of J^''N while OAl and OA2 include those of ]^^0. We note from 
table|5]that the minimum values obtained in the alternative fitting procedure are even 
lower than than those shown in table [3] with the set OAl having the lowest value 
among all the sets. We have used the parameter sets HPA2 (that corresponds to min- 
imum;^f^ in table [3]), NAl and OAl of table |5] in our calculations on the hypernuclear 
properties presented in the next section. 

4. Results and Discussions 

4.1. A single particle fields 

The existence of the ground states of A hypernuclei is related to the single particle 
field (Va) for the A hyperon. In the SHE model the AA^ force contributes to the field 
Va in proportional to various powers of p^- It would be instructive to compare this 
quantity calculated within the Hartree-Fock theory (V^'') with those obtained with a 
phenomenological parameterization of (V^'') e.g. of Ref. 129] for hypernuclei with 
masses varying in a broad range. In Eig. [U we show [calculated using Eq. ( fT3T l 
with the parameter set HPA2], for hypernuclei O, ^°Ca, ^''Y, ]^'^La and ^o^Pb. We also 
show in this figure the V^'' obtained within the parameterization of Ref. [29] for O, Ca 
and Pb systems - the corresponding potentials are actually taken from this reference. 
We notice that the central depths of potentials V^'^ are around 29 MeV irrespective of 
the nuclear mass. In comparison, V^* are somewhat less attractive. This result is in 
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2 4 6 8 10 
r (fm) 



Figure 1 : [color online] The solid lines show the coordinate space SHF mean field potential Va in hypernuclei 
]^0, ^"Ca, fY, ^ha, and ^ospb as marked near each hne. The corresponding fields generated in an empirical 
Woods-Saxon parameterization (that are taken taken from Ref. (2Sj| ) are also shown for O, Ca and Pb nuclei 
by dashed lines near the lines showing the corresponding Hartree-Fock Va- 

contrast to that shown in Ref. ifsoll where the V^'' was found to be shallower than V^' . 
Use of a different AN force could be one reason for this difference. 

For the Pb nucleus, the extends to the larger values of r in comparison to V^''. 
Similar trend was seen in Ref. ll39ll . This result has been attributed to the presence 
of the effective mass m* in the HF equation which mocks up for some finite range ef- 
fects in such approaches. This results is also in agreement with the observations made 
in Ref. 1E9I1 where it was shown that the presence of the terms in the A potentials 
[Eqs. (13) and (14)] leads to an interaction with an increased half value radius as com- 
pared to that of the simple Woods-Saxon potential for heavier nuclei. However, in 
lighter systems the dependent term of the potential seems to reduce the effective 
radius of the potential. 

4.2. Binding energies of A single-particle states 

In Fig. 2, we show a comparison of the calculated and the experimental A binding 
energies (Ba) of Is, Ip, Id, I f and Ig shells of a number of hypernuclei with masses 
ranging from 8 to 208. In these calculations, we have used SLy4 for the A^A^ interaction 
and the parameter sets HPA2, NAl and OAl for the AA^ interaction. We note that 
there is an overall agreement between our calculations performed with HPA2 set and 
the experimental data in both the low mass and the high mass regions. However, sets 
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Figure 2: [color online] The comparison of experimental and calculated binding energies of Is, \p, \d. If, 
and Yg A states as a function of A"-'' where A is the baryon number of the hypernucleus. In each case the 
solid, dashed and dashed-dotted lines coiTespond to the results obtained with A7V interaction parameter sets 
of HPA2, NAl and OAl, respectively. 

OAl and NAl tend to overbind the Is and also to some extent the \p A orbitals for 
lighter nuclei. This effect is relatively more stark with the set NAl. However, since 
the set HPA2 has been obtained by including the BE of the light nuclei as well in the 
fitting procedure, it is not surprising that it provides a better description of the Z?a of 
the \s and \p orbitals of the lighter systems. 

For heavier nuclei all the three sets describe the Ba equally well. In a marked con- 
trast to the SHF results of Ref . [ 39] , we observe the underbinding only for the 1 s single 
particle orbital of ^'^*'Pb hypernucleus. Moreover, even for this case the magnitude of 
the underbinding is much smaller in our calculation as compared that of Ref. [39]. In 
the SHF calculations presented in Ref. 1,46] , underestimations of varying magnitudes 
were observed for the single particle energies of both s and p orbitals of the heavier 
systems with most of the AA^ interactions used. In the density dependent hadron field 
theory calculations of Ref. |5], the single particle energies were underpredicted by fac- 
tors of up to 2.5 for hypernuclei with masses below ^**Si. This was attributed to the 
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Table 6: The A effective mass (calculated at saturation density po) and rms radii of ground state A orbitals 
of various hypernuclei. 



SET 


my Ota 


rAilLi) 






rAifSi) 


rAQ'Ca) 


rA(^'V) 


rA^^'La) 








(fm) 


(fm) 


(fm) 


(fm) 


(fm) 


(fm) 


(fm) 


(fm) 


HPAl 


0.88 


3.16 


2.96 


2.98 


3.09 


3.27 


3.42 


4.10 


4.75 


HPA2 


0.85 


3.17 


2.98 


3.00 


3.12 


3.31 


3.46 


4.17 


4.81 


NAl 


0.88 


2.94 


2.88 


2.92 


3.11 


3.31 


3.46 


4.14 


4.81 


OAl 


0.99 


2.94 


2.78 


2.82 


2.93 


3.11 


3.27 


3.87 


4.54 


HPA3 


0.86 


3.19 


2.98 


2.99 


3.09 


3.27 


3.43 


4.11 


4.75 


HPA4 


0.87 


3.19 


2.98 


3.00 


3.10 


3.28 


3.43 


4.12 


4.75 


NA2 


0.88 


2.95 


2.87 


2.92 


3.09 


3.28 


3.44 


4.09 


4.76 


OA2 


1.00 


2.93 


2.74 


2.80 


2.90 


3.07 


3.23 


3.82 


4.49 



enhancement of surface effects in the lighter nuclei that are not accounted for properly 
in their calculations. 

Looking at the fact that the parameter sets HPA2 provides a good overall descrip- 
tion of the single particle spectra of the A hypernuclei which in some ways improves 
upon previous SHF calculations performed with different parameterizations, the pre- 
dictions of our model can be used with relatively more confidence in calculations of 
the hypernuclear production cross sections via e.g. (7t^,K'^) fl^], or (y,K'^) lITsifTill 
reactions for those cases where experimental information may not be available about 
the binding energies of the A states. 

In table|6] we show the effective A hyperon masses and the root mean square (rms) 
radii (ta) of the ground state A orbits of various hypernuclei calculated at the nuclear 
saturation density («o) with various AN parameterizations obtained by us. This pro- 
vides additional constraints on different parameter sets. First of all we note that m*JmA 
obtained with interactions having y value of (1/3) (last 4 rows of this table) are slightly 
different from those involving 7=1. The effective mass is largest with the parameter 
sets OAl and OA2. Furthermore, these sets give the lowest ta for various hypernuclei. 
It is also of interest to note that ta calculated with our parameter sets are systematically 
larger than those obtained with Skyrme AA^ interactions of Ref. [46]. It should however 
be mentioned that some uncertainty in the rms radii originates from the choice of the 
NN potential which was also different in Ref. 1.461 . 

4.3. Total binding energy per baryon of hypernuclei 

Although the hyperon contribution to the total hypernuclear energy is relatively 
small, a systematic study of the total binding energies per baryon {BE /A) of hypernu- 
clei could be useful because it sheds light on the pattern of the stability of hypernuclei 
across the periodic table. Looking at the experimental data one notices that the BE/A 
of the medium mass hypernuclei is larger than those of the lighter and the heavier ones. 
This indicates that hypernuclei in this region may be more stable. In the SHF frame- 
work, the total binding energy of a hypernucleus of baryon number A can be obtained 
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Table 7: Binding energy per baryon number (BE/A) of vaiious known single lambda hypemuclei and con'e- 
sponding values calculated by RMF method |79] are also listed for comparison. 



Hypernuclei BE/A BE/A (RMF) Hypernuclei BE/A BE/A (RMF) 
(MeV) (MeV) (MeV) (MeV) 



\l 


7.115 


- 


fNi 


8.773 


8.839 


5.906 


- 


8.773 


8.844 


8.631 


- 


ANi 


8.768 


8.853 


aN 


7.076 


- 


aNI 


8.764 


8.865 


JfN 


8.523 


- 




8.727 


8.778 


^"Ne 


7.386 


7.635 


a'Kt 


8.695 


8.796 


i'Ne 


7.579 


7.714 




8.745 


7.680 




7.624 


7.723 


^^Rb 


8.719 


7.587 


f Al 

33 o 


8.184 


- 


88 c 

89 c 
A 


8.754 


7.657 


8.666 


- 


8.734 


7.564 


8.579 


8.637 


A ^ 


8.739 


7.510 




8.579 


8.704 




8.740 


7.545 


3SAr 


8.585 


8.597 




8.731 


7.455 




8.634 


8.728 


^2Mo 

f'Mo 


8.675 


7.378 


toTi 


8.863 


8.855 


8.678 


7.291 


8.749 


8.880 


-Sn 


8.471 


6.925 


8.823 


8.843 




8.479 


6.961 




8.747 


8.904 


i-Sn 


8.487 


6.996 


52y 


8.751 


8.909 


-Sn 


8.495 


7.032 


8.725 


8.774 




8.502 


7.030 


fNi 


8.715 


8.897 


i-Sn 


8.509 


7.028 


8.666 


8.856 


-Sn 


8.516 


7.027 


8.705 


8.847 


rsn 


8.523 


7.026 




8.738 


8.841 


8.530 


7.026 




8.763 


8.838 




8.538 


7.026 



from Eq. (22). 

In Tables (|7]i and (O we present the results of our SHE calculations (done with the 
HPA2 force) for BE/A for 73 hypernuclei with masses in the range of 9 to 21 1. The 
results obtained with sets NAl and OAl are almost the same. We note that for the 
lighter systems (A < 28), the stability of a hypernucleus depends on that of the nucleus 
where a neutron is replaced by the hyperon. For example, ^Be is more stable than ^B, 
and is more stable than '^N. For hypernuclei with baryon number in the range of 
30-95 the BE /A is around 8.7 MeV and it decreases gradually with further increase in 
the mass number. For comparison purpose the results of a relativistic mean field (RMF) 
model taken from Ref. 1,791 are also presented wherever available. RMF results for the 
BE /A are given only in this reference. We note that that for nuclei with baryon number 
in excess of 87 the RMF model underbinds the hypernuclei. However, the difference 
seen between SHE and RMF energies could possibly result from the differences in the 
nuclear energies predicted by the two models. Therefore, one should be careful in 
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Table 8: Same as table [7] 



Hypernuclei 


BE/A 


BE/A (RMF) 


Hypernuclei 


BE/A 


BE/A (RMF) 




(MeV) 


(MeV) 




(MeV) 


(MeV) 




8.545 


7.026 




8.204 


6.874 




8.551 


7.027 




8.337 


6.823 




8.554 


7.028 


I42p, 


8.339 


6.855 




8.552 


7.030 




8.321 


6.800 




8.415 


7.066 




8.325 


6.834 




8.391 


7.074 


rsm 


8.280 


6.753 




8.392 


6.949 




8.289 


6.790 




8.381 


6.968 


7.886 


6.726 


A 


8.375 


6.886 


2iOBi 


7.854 


6.711 


8.369 


6.911 


7.826 


6.689 


f'La 


8.362 


6.893 


7.833 


6.696 




8.352 


6.845 









interpreting the comparison of the results produced by these two mean field models. 



4.4. Properties of a neutron star 

In order to calculate neutron-star properties, it is necessary to have an Equation 
of State (EOS) linking pressure to the total energy density of the dense matter. At 
densities closer to the nuclear saturation density («o = 0.1 6/m"-'), the matter is mostly 
composed of neutrons, protons and leptons (electrons and muons) in fi equilibrium. As 
density increases, new hadronic degrees of freedom may appear Hyperons are one of 
them as the equilibrium conditions in neutron stars make the formation of hyperons 
energetically favorable. The role of hyperons on the neutron sta r prope rties has been 
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studied by several authors (see, e.g. Refs 

In this section, we employ the best fit AA^ interactions obtained by us to discuss the 
implications of hyperons on the EOS and the structure of neutron stars. Unlike several 
other authors, we have included only the A hyperons into our calculations as this study 
is restricted to testing the interactions obtained in this work. In that sense our work may 
appear to be less complete in comparison to many previous studies where include other 
hyperons were included as well. Nevertheless, in all likelihood the SA^ interaction is 
repulsive because no stable Z hypernucleus other than that of mass 4, is known to exist. 
Therefore, S appears at much higher densities ll87ll as compared to A. The H hyperon, 
on the other hand, could appear at densities comparable to those of A. However, there 
is considerable amount of uncertainty about the strength of "EN interaction as no bound 
E hypernucleus has been detected so far. The threshold of the appearance of the H 
hyperons is pushed to higher densities with increasing 'B.N potential. 

In our calculations of the EOS, we have closely followed the methods reported in 
Refs. ilS 1471 14911 . To our energy density functional [Eq. (1)] we have also added a 
contribution corresponding to the AA interaction as at higher densities the results are 



sensitive to this term, for which the energy density functional is taken from Ref. 147 1 
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Figure 3: [color online] The equation of state (pressure as a function of baryon number density) obtained 
by using AN interactions HPA2 (dashed line), HP A3 (dashed-dotted line), NAl (dashed-double dotted line) 
and OAl (dotted line). 



with parameters corresponding to their set SLL2. The energy density includes the rest 
energies of the matter constituents. The leptonic contribution to the energy density is 
calculated as discussed in Refs. |47, 4^ (see also [91]). For this sector the energy 



densities corresponding to electrons and muons are written as 

1 r*/ I 

S, = — I k^^k'^+m]dk, (25) 



where I - e (electron) or /i (muon) and A:^ is the corresponding Fermi momentum. 
For a given baryon density the values of the Fermi momenta for neutrons, protons, 
A, electron and muon can be obtained by requiring that neutron star matter is in p 
equilibrium. This leads to the following equations for the equality of the chemical 
potentials (represented by n in the following) 

H„ -flp = jUe, IJf, = fie 

fi„ + m„ = jUa + niA, (26) 
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Figure 4: [color online] Particle fractions obtained with using AA' interactions HPA2 (a), HPA3 (b), NAl 
(c) and OAl (d). 



where chemical potentials are defined as 



5e 



(27) 



where & is total energy density and rij the particle number density. 

The total baryon number density is ni, - + Up and the charge neutrality requires 
Up - He + tift, where «f and «^ are the number densities of electrons and muons, respec- 
tively. These equations combined with Eq. |26]give the particle fraction Yj = ^. The 
EOS is defined by the expressions 



(28) 



where p{nb) is the mass density of the matter 
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Figure 5: [color online] Mass radius relation for neutron stars obtained with EOS shown in Fig. 3 with 
interactions HPA2, HPA3, NAl and OAl. 

In Fig. 3 we show the equation of state calculated with AN interactions HPA2, 
HPA3, NAl, and OAl. It is to be noted that in each case the inclusion of hyperons 
makes the EOS much softer with respect to that of the pure nucleonic case. Since hy- 
perons can be accommodated in the lower momentum states, their kinetic energies are 
decreased which leads to the softening of the EOS. We see that the degree of softness 
of the EOS obtained with HPA2 and NAl interactions are almost identical. However, 
with HPA3 the softness is comparatively smaller and with OAl the softening of the 
EOS is relatively the lowest. 

In Fig. 4, we show the particle fraction Yj as a function of the total baryonic density. 
It is clear from this figure that neutrons are by far the most dominant object of the 
matter One further notices that while in panels (a), (b) and (c) A hyperon appears at 
about twice the nuclear saturation density, in panel (d) they start only after about 2.5«o. 
They are also less numerous in this case as compared to those in other 3 panels. This 
is the consequence of the fact that with the interaction OAl the softness introduced to 
the nucleon only EOS is the least. The rise of hyperon fraction, after their formation, 
is quickest with HPA3 and NAl interactions as compared to that seen with with other 
interactions. 

Further consequences of the different nature of the EOS softening with the intro- 
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duction of the A hyperon with different AA^ interactions can be seen on the their effect 
on the neutron star masses. To obtain the relation between neutron star mass and its 
radius, we have solved the Tolmann-Oppenheimer-Volkoff equation 

dP GS(r)M(r) r ^ P(r)_^r ^ AnPPjr) -, 



dr r^c^ L- £(r)JL' " M(r)c^ 

2GM(r)n 



and 



dM{r) _ Anr^&ir) 
dr 



(30) 



where P, S, and M are the pressure, energy density and gravitational mass of the 
neutron star, respectively. These quantities depend on the distance r from the center. 
Eqs.|29]and[30]can be solved from the knowledge of the initial value of the pressure P 
at the center (r =0) and using M(0) - and the relationship of P and & (EOS). 

The predicted masses of the neutron star [measured in solar masses Mq] as a func- 
tion of radius (in kilometer[km]) are shown in Fig. 5 using EOS of Fig. 3. We note that 
with EOS obtained with interactions HPA2, HPA3, and NAl maximum neutron star 
mass is similar (around 1.5 Mq). In Fig. 5, Mq is depicted as Mj^- However, with the 
EOS corresponding to the interaction OAl the maximum mass is about 1.75 Mq. This 
result can be understood from the fact that a stiffer EOS leads to a larger neutron star 
maximum mass. It may be mentioned here that if AA interaction is switched off the 
maximum neutron star mass fails to reach to even its "canonical" value of 1.4 Mq. 



5. Summary and conclusions 

In summary, we have used an extended Skyrme Hartree-Fock model to describe the 
properties of the A hypernuclei. New parameterizations for the Skyrme type AN force 
have been obtained by fitting to the experimental binding energies of about 20 A hyper- 
nuclear orbitals with the baryon number ranging from 8 to 208 (the best fit parameter 
set is termed at HPA2). We have also performed the fittings by excluding the binding 
energies of nuclei with masses below 16. In this case two sets of parameters were ob- 
tained: in the first one (termed as NAl) the binding energies of the ]^^0 hypernucleus 
were excluded but those of the J^^N hypernucleus measure recently in a (e, e'K"^) experi- 
ment were include in the fitting procedure, while in the second one (termed as OAl) the 
reverse was done. The fitting method uses an elegant minimization method based 
on the simulated annealing method. 

These three sets of the best fit parameters were used to calculate the binding ener- 
gies of 1 s, Ip, Id, If and Ig shells of a number of hypernuclei. We find that calcula- 
tions performed with the set HPA2 provide the best agreement with the experimental 
binding energies in the entire range of mass values and for all the orbitals. The set 
OAl produces binding energies, which overestimate somewhat the corresponding ex- 
perimental data of the 1 s and to a lesser extent of the 1 p orbitals for the lighter nuclei. 
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This overestimation is larger in case of the set NAl. However, all the three sets pro- 
duce equally good agreement with the experimental data for heavier nuclei for all the 
orbitals. 

Except for one case [1 5 orbital of ^"'^Pb nucleus], our calculations do not underbind 
the heavier systems. This is a marked improvement over the similar previous Hartree- 
Fock studies of A hypernuclei where binding energies of orbitals of several heavier 
nuclei were underpredicted. Furthermore, the root mean square radii of heavier systems 
are predicted to be larger than those obtained in previous SHF calculations done with 
different AA^ forces. 

We made a systematic study of the mass dependence of the total binding energy 
per baryon of 73 A hypernuclei spanning the entire range of the periodic table. It is 
observed that hypernuclei with masses in the range of 30-95 are more stable than those 
lying in other regions. For lighter systems some of the hypernuclei are more stable than 
their immediate neighbors on both sides. 

We have also tested our best fit AA^ interactions to investigate the role of hyperons 
in the neutron star sector In these studies, we used the same A^A^ effective interaction 
(S Ly4) together with HPA2, HPA3, NAl and OAl effective AN interactions. Inclusion 
of the AA interaction was also found to be necessary in these studies - we took the 
corresponding effective Hamiltonian from the work of Lanskoy ll94tl . It is noted that 
inclusion of the A hyperon makes the neutron star equation of state softer However, 
while the sets HPA2, HP A3 and NAl lead the softness of the similar magnitude, that 
obtained with the set OAl is lesser Furthermore, with this set the A hyperon appear at 
a relatively higher density and are relatively less numerous as compared to other three 
sets. 

The maximum neutron star mass obtained with sets HPA2, HPA3 and NAl is about 
1.5Mo, while that with the set OAl is I.TSMq. This is the direct consequence of the 
relatively stiffer EOS in this case. We remark that these results obtained with our best 
fit AA^ effective interactions reproduce quantitatively all the features that have been 
observed in several different models such as non-relativistic Brueckner-Hartree-Fock 
and relativistic mean field calculations. Nevertheless, an explanation of some of the 
recent neutron star observations 19211 such as pulsar PSR J1614-2230 with a mass of 
I.97+O.O4M0 is still beyond the scope of the present Skyrme type of model unless the 
EOS is even stiffer than what has been achieved with the interaction OAl. 
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